%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Fonction d'assemblage du vecteur effort généralisé par élément
%
% Entr�es :
%
% sous-maillage : vecteur de coordonnées des noeuds du sous-maillage
% uh : [vecteur 4x1] extraction de la solution éléments finis
% epsilon : =1 sur le premier �l�ment du patch,
%           =-1 sur le deuxi�me �l�ment du patch,
%
% Sorties :
%
% ie : [matrice 4x4] matrice ie

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [ Fe ] = Assemblage_F_element_fluxfree( maillage , uh , epsilon , poutre,sousPoutre,donnees )


nse = size(maillage,1)-1 ;                % nombre de sous-éléments
K = zeros(2*nse+2,4) ;


for j=1:nse % attention : i a �t� chang� en j !!

%     [a,b,c,d]=coeffs_abcd(i,ne,epsilon) ;
%     ie(:,2*(i-1)+1:2*(i+1)) = ie(:,2*(i-1)+1:2*(i+1)) + Ie(a,b,c,d,epsilon,poutre,sousPoutre) ;

[a,b]=coeffs_ab(j,nse);

% D�finition de la "matrice de rigidit�" associ�e au calcul du vecteur
% effort : Kf=Kf1+Kf2

K(2*(j-1)+1:2*(j+1),:) = K(2*(j-1)+1:2*(j+1),:) + Kf(poutre,sousPoutre,epsilon,a,b,nse,(donnees.nombreElements)) ;




end

Fe=(K*uh) ;


end

